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1 Introduction 
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We study the U(2) lattice gauge theory in 
the pure gauge sector using the simplest 
action, with determinant and fundamental 
terms, having the naive continuum limit of 
SU(2)xU(l). We determine part of the phase 
diagram of the model and find a first-order 
critical line which goes through the U(l) crit- 
ical point. We show how to deduce both 
the order parameter of the first-order transi- 
tion and the U(2) renormalization group fiow 
from the lattice potential in the determinant 
and fundamental representations. We give ev- 
idence that a Monte-Carlo simulation of the 
model is indeed consistent with the symmet- 
ric SU(2)xU(l) continuum limit in the weak 
coupling pertubative regime. 
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A description of the weak and electromagnetic 
interactions by a local gauge theory based on 
the group U(2) is consistent with the quan- 
tum numbers of the observed elementary par- 
ticles. Indeed the global structure of the stan- 
dard model can be inferred from the observed 
patterns of the irreducible representations of 
matter fields which obey the quantization rule: 
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where T and Y are the weak SU(2) isospin and 
weak U(l) hypercharge respectively and Q is 
the electric charge. It has long been empha- 
sized [1] that the global structure of groups is 
relevant to particle physics. 

On the lattice, the formulation of the 
U(2) gauge theory differs from that of the 
SU(2)xU(l) gauge theory already in the pure 
gauge sector. Indeed U(2) is isomorphic to the 
quotient group 
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SU(2) X U(l) 
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which means that each element \J of the group 
U(2) can be represented as a pair (e"^,T/) G 
U(l)xSU(2), with the identification (e''^,y) = 



(- 



-V). Wilson's lattice gauge invariant 



formulation [2] depends explicitly on the topo- 
logical structure of the gauge group because 
the degrees of freedom belong to the group 
manifold, not its Lie algebra. It is very often 
stated that compactness of the group mani- 
fold is a lattice artifact because it could in- 
duce spurious effects in the approach to the 
continuum limit. In fact a global formulation 
of gauge theories also exists in the continuum 



1 



in terms of loop variables [3]. Such a non-local 
description of gauge theories contains all their 
kinematical properties, though it is not obvi- 
ous how to write down a classical action for 
loops. We shall consider the possibility that 
the global structure of gauge groups might 
play a role in the non-pertubative properties 
of gauge theories. But we must impose that 
any admissible lattice formulation does repro- 
duce the results of the continuum non compact 
formulation in the pertubative regime. 

There is a large degree of arbitrariness in 
the choice of action on the lattice. In particu- 
lar the most general U(2) lattice action which 
can be built out of the elementary plaquette 
P reads 

S[Up] = J2f^^ Tr^ Up + li-c- 

r 

= ^/3n,ie'"^-Tr,-T/p+h.c.(3) 

"J 

where the sum runs over all irreducible rep- 
resentations (irreps) r of U(2), which are the 
direct products of the U(l) irreps (labelled by 
an integer n ^ Z) and SU(2) irreps (labelled 
by a non-negative integer or half-integer j), 
subject to the constraint (— = 1^ which 
is a consequence of Eq. (1). 

It is generally believed that the choice of ac- 
tion is largely irrelevant once the gauge group 
is fixed, because it amounts to a modifica- 
tion of the regularization scheme which should 
not affect the continuum limit as long as weak 
coupling pertubation theory is applicable. In 
particular the U(2) lattice gauge theory is ex- 
pected to reproduce the SU(2)xU(l) asymp- 
totic scaling properties. 

If we parametrize an SU(2) group element V 
as a 2 X 2 unitary unimodular matrix, the Wil- 
son action for the SU(2)xU(l) lattice gauge 
theory without matter fields reads: 

S[U] = J2hcosct,p + ^Ti(Vp + V^) (4) 
p 

Because of the decoupling between the U(l) 
and SU(2) gauge fields, its properties can be 
directly inferred from those of the U(l) and 
SU(2) lattice gauge theories with Wilson ac- 
tion, which are well-know [4, 5], and are sum- 
marized in Fig. 1. There is a U(l) first-order^ 
critical line at yf ~ 1.01 and the renormaliza- 
tion group trajectories are the straight lines 
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Figure 1: Critical curve and renormalization 
group trajectories of SU(2)xU(l) lattice gauge 
theory with Wilson action (solid lines) and ex- 
pected results for U(2) lattice gauge theory 
with extended action (dashed lines). 



parallel to the &2 axis in the U(l) cold phase 
hi > hf. 

They are labelled by a constant effective 
U(l) charge e(&i). In the compact formulation 
there is a natural dual description of mass- 
less free photons in terms of monopoles loops 
which explains the renormalization of the bare 
charge hi in the cold phase and the "con- 
finement" transition through monopole con- 
densation [11]. Along each renormalization 
group trajectory the physics is that of a con- 
fined SU(2) gauge theory which fixes the scale 
through the string tension. In the U(l) hot 
phase hi < hf there is no continuum limit. 

The simplest actions for the U(2) lattice 
gauge theory which have a naive continuum 
limit consistent with the two bare couplings 
of the SU(2)xU(l) gauge theory are obtained 
from Eq. (3) by retaining the first two terms 
consistent with the condition (— 1)"+2J = 1^ 
namely n = 2, j = and, n = 0, j = 1 or 
n = l,j = 1/2: 



^The transition is first-order at least for periodic 
boundary conditions, i.e. wlien tlie lattice is a torus 

J-<4: 
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S[U] = hi COS 2<t)p 

p 

+ yTrl/pTrl/; (5) 

S[U] = ^ &icos2^Ap 
p 

+ ^-^i:i{e'-^^Vp + e-'-^^vl) (6) 

Eq. (5) is notliing but a lattice action for 
the SO(3)xU(l)/Z(2) gauge theory. The de- 
coupling between the U(l) and S0(3) gauge 
fields on the lattice makes its asymptotic scal- 
ing properties in the continuum similar to 
SU(2)xU(l). However the S0(3) gauge the- 
ory does not exhibit confinement since there 
is a first-order critical point at hij ~ 2.5 [6] 
and the S0(3) lattice potential is screened at 
large distances. We have here a noteworthy 
example of the dependence of the non pertu- 
bative properties of a gauge theory upon the 
topological structure of its gauge group. 

On the other hand action (6) introduces 
a coupling on the lattice between the U(l) 
and SU(2) gauge fields. It is usually expected 
that this coupling disappears in the continuum 
limit. The naive continuum limit of this action 
is, up to an additive constant: 

X, (l-ll') 

+ '-fTrW^^Ax) (7) 

where a is the lattice spacing, Ffi^(x) and 
Wi,^(x) the U(l) and SU(2) field strength ten- 
sors respectively. In this limit the renormal- 
ized group trajectories in the (hi, 62) plane are 
the straight lines with equations: 



They are sketched in Fig. 1. The analogy with 
SU(2)xU(l) would imply that these lines are 
labelled by an effective coupling e(gi) to a 
massless U(l) gauge field. Along these lines 
the physics should again be the same as for a 
confined SU(2) gauge theory. 

The U(2) lattice gauge theory has been 
studied [7, 8] with the Wilson action, which 
corresponds to letting hi = in Eq. (6). It is 
therefore not clear how to extract the contin- 
uum limit from these studies which were done 



in the early days of Monte-Carlo simulations 
of lattice gauge theories. However a convinc- 
ing signal for a first-order transition was found 
for &2 ~ 3.3 and some evidence was found 
for a decoupling between the U(l) and SU(2) 
gauge fields. This decoupling had also been 
established for the Wilson action of U(N) lat- 
tice gauge theory in the large N limit [9]. It 
was clearly recognized [7] that a determinant 
piece should be added to the Wilson action of 
U(2) lattice gauge theory^ but the computing 
power was not available at the time. To our 
best knowledge this problem has not been ad- 
dressed since then. 

Our purpose is to report on a study of U(2) 
lattice gauge theory with extended action (6). 
In the next section we make a prerequisite step 
which is the determination of the phase dia- 
gram of such a lattice model. In this work we 
describe only the phase diagram in the first 
quadrant (&i > 0, &2 > 0) of the parame- 
ter plane. We find a first-order critical line 
which joins, as anticipated in Fig. 1, the points 
A (1.0, 0.) and 5(0., 3.3). We then confirm 
[7, 8, 9] that the determinant of the Wilson 
loops can serve to characterize the transition. 
The determinant probes the central U(l) sub- 
group of U(2) and we can extract its lattice 
potential. We find that one phase is con- 
fined with a non-zero string tension and the 
other phase is deconfined with an effective 
U(l) charge. We determine the contour levels 
of these effective charges and find that, asymp- 
totically, they are straigth lines, in agreement 
with the naive continuum analysis. We also 
study the lattice potential in the fundamental 
representation of U(2). In the deconfined U(l) 
phase we find a non-zero string tension which 
it is natural to associate to the SU(2) part of 
U(2). We study the renormalization fiow of 
the bare coupling of the SU(2) subgroup by 
holding fixed the string tension and find that 
it is consistent with the contour levels of the 
effective U(l) charge and with universality. 

In our conclusion we make some remarks 
about the limitations of this work and further 
possible extensions. Let us mention here that 
the phase diagram has also a rich structure 
in the quadrant (&i < 0, &2 > 0). The well- 

^ Samuel has also made an analytic study of the 
large N limit of U(N) lattice gauge theory with funda- 
mental and ajoint action terms [10] , which corresponds 
to retaining the terms n = 1, j = 1/2 and ri = 0, j = 1 
in Eq. (3). 
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known hysteresis phenomenon at first-order 
transitions makes difficult a reasonably pre- 
cise determination of this phase diagram. It 
remains to unravel whether the additional fea- 
tures are relevant to the continuum limit of the 
U(2) lattice gauge theory and how the phase 
structure is related to the monopole content of 
the model [12]. 



2 Phase diagram 

In order to determine the phase diagram of 
the U (2) lattice gauge theory with action (6) 
a standard Monte-Carlo simulation has been 
done on an S'* lattice with periodic boundary 
conditions using a local heatbath algorithm. 
There is, perhaps, one technical point worth 
mentioning in the implementation of the heat- 
bath algorithm. Because action (6) contains 
both determinant and trace terms, it proves 
more convenient to use a parametrization of 
the U(2) group which incorporates its quotient 
structure: 



U 



Ul + tU'j 

'(-"3 + iui) 



"3 + 



t U4 



(9) 



with u-^ + U2 + + = 1 and < 9 < 2tt. 

Indeed with this parametrization action (6) 
depends linearly upon cos6', sm9 and the Ui 
variables. The U(2) Haar measure is still a di- 
rect product of the form dlJ ~ dOd'^ubiu'^ — 1). 
Therefore we can generate a U(2) group ele- 
ment by using standard or optimized heatbath 
algorithms for U(l) and SU(2). 

This study has been restricted to the first 
quadrant (&i > , &2 > 0) in the parameter 
plane. Data have been taken for several hun- 
dreds of points with 2000 iterations for each 
point after a cold start. Each run has been 
analyzed by the histogram technique with the 
PAW package [22]. In order to determine the 
phase diagram we have included in the anal- 
ysis the measurements of the trace and deter- 
minant of the elementary plaquette. Average 
values and standard deviations are obtained 
by a gaussian fit to the histograms. Unther- 
malized data are spotted by a bad of the 
fit. Typically a few hundred of iterations are 
discarded for thermalization. The number is 
adjusted from the . First-order transitions 
are located by looking for histograms with a 
two peak structure fitted by two gaussians. 
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Figure 2: Typical histograms in the critical 
region. 



Typical two-peak histograms and their fit are 
shown in Fig. 2 for the trace and determinant. 

We have a very clear double-peak signal for 
a first-order transition in the {hi, 62) plane. 
Critical couplings are defined by fits with two 
gaussians of equal weight. The critical cou- 
plings coincide for the trace and the determi- 
nant. Their values are displayed in Table 1. 
The first-order critical curve is shown in Fig. 3. 
It joins the U(l) critical point at (1.0,0.) to 
the U(2) critical point at (0.,3.3) [7]. 

A by-product of analyzing first-order tran- 
sitions with histograms is the direct measure- 
ment of the discontinuities of observables at 
the transitions. The gaps in the plaquette en- 
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0.048 


0.99 


0.5 
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0.045 


1.01 


0.0 




0.0 


0.045 



Table 1: Critical couplings of U(2) 8^* lattice 
gauge theory after a cold (c) and hot (h) start. 
AT and AD denote the jumps in the trace and 
determinant of the plaquette across the critical 
line. 




Figure 3: First-order critical curve of U(2) 8^* 
lattice gauge theory. Solid line corresponds to 
cold start, dashed line to hot start. Lines in- 
terpolate linearly between data points denoted 
by markers. 
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Figure 4: Gaps in the plaquette energy along 
the first-order critical curve for U(2) 8^* lattice 
gauge theory. 



ergies can be read from Table 1 which con- 
tains the values of the jumps in the trace and 
determinant of the plaquette along the criti- 
cal curves in the cold and hot phases. These 
gaps are also plotted in Fig. 4. We see clearly 
two regimes in the gaps if we follow the criti- 
cal curve labelled by the coupling hi. Starting 
from the U(l) critical point at hi = 1.0 the 
gap in the determinant stays rather constant 
down to hi ~ 0.7 where it departs from the 
U(l) value. 

There is a well-known problem that we must 
worry about when studying first-order transi- 
tions which is the hysteresis phenomenon. The 
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Figure 5: Typical hysteresis cycle in the criti- 
cal region. 



autocorrelation time for local Monte-Carlo al- 
gorithms is very long at a first order transi- 
tion because any local algorithm requires the 
system to pass through the minimum sepa- 
rating the two peaks. This tunneling time 
increases exponentially with the gap between 
the peaks. The aforementionned widening of 
the gap makes the critical slowing down even 
worse for U(2) than for U(l). In order to assess 
the problem more quantitatively we have re- 
peated the Monte-Carlo runs with a hot start 
for all points in the critical region. Everything 
else is unchanged: local heatbath algorithm, 
2000 iterations, histogram analysis. The cor- 
responding critical couplings are shown in Ta- 
ble 1 and plotted in Fig. 3. We observe that 
the hysteresis A& = — correlates well 
with the width of the gaps. A typical hystere- 
sis cycle is shown in Fig. 5. 

These hysteresis effects make a precise de- 
termination of the critical curve practically 
impossible with a local algorithm. The stan- 
dard prescription to locate the phase transi- 
tions is to make runs from mixed starts which 
corresponds roughly to take the average of the 
two values listed in Table 1. The hysteresis 
induces systematic errors in the values of the 
gaps of Table 1 which are much larger than 
the statistical errors. We estimate these sys- 



temactic errors to be of the order of 0.01. In 
the sequel our focus will be rather on the bulk 
properties of the phases which do not depend 
upon the precise location of the critical curve. 



3 Order parameter 

The critical branch found in the previous sec- 
tion goes through the U(l) critical point at 
hi = 1.0, &2 = 0. Since only the determinant 
term in the action plays a dynamical role in 
the vicinity of this point, it is natural to expect 
that the order parameter of the phase transi- 
tion of the U(l) subgroup probed by the de- 
terminant is also the order parameter of the 
U(2) phase transition along this branch. The 
order parameter of the U(l) phase transition is 
the string tension [4] which is zero in the cold 
phase and non-zero in the hot phase. There- 
fore we need to extract the string tension pro- 
duced by the central U(l) subgroup of U(2) 
from the determinant of Wilson loops. 

We use a method [13] which is very effec- 
tive on intermediate size lattices and would be 
quite useful also on large lattices. In Monte- 
Carlo measurements of the string tension, one 
needs operators on the lattice which are usu- 
ally built either from planar rectangular loops, 
or from thermal lines'*. In both cases an ex- 
trapolation in one, say the time, direction is 
necessary to make contact with the physical 
observables. This extrapolation introduces an 
important source of systematic errors. There 
is a convenient way to circumvent this prob- 
lem. It consists of building lattice operators 
which contain only one length scale. This cri- 
terion will be met if one uses loops inscribed 
on an hypercube. It is easy to form, with 
one-scale non-planar loops, ratios from which 
ultraviolet divergences can be removed by a 
charge renormalization only: 



< det W(C, L) > 
< detW(L,L) >P/^ 



(10) 



where /3 denotes, collectively, the coupling pa- 
rameters on the lattice, det W(C, L) is the 
determinant of the U(2) phase factor along 
the loop C inscribed on an hypercube of side 
La and AeiW{L, L) is the determinant along 
the square loop of side La to the power P /A 



^For a recent review about string tension see [14] 
and references about standard algorithms therein. 
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Loop 


Description 


SiC) 
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1, 2, 3,-1,-2,-3 


1.930 


0.394 


C2 


1, 2, 3, 4,-1,-2,-3,-4 


3.014 


0.714 


C3 


1, 2, 3, 4,-2,-1,-4,-3 


2.899 


0.525 



where PL is the perimeter of C in lattice spac- 
ing units a. Then all perimeter and corner 
singularities are removed [15]. Wilson's area 
law gives the leading contribution to the ratio 
Dc{L,j3) and we can write, from dimensional 
analysis alone: 

Dc{L,l3) = exp{-KiiSc -P/i)L''a' 
-Bc + ---} (11) 

where A'l is the string tension in the determi- 
nant representation, Sc is the area of the min- 
imal surface enclosed by the loop C inscribed 
on the unit hypercube. Be is a constant re- 
lated to the effective charge of the Coulomb 
part of the potential and the discarded terms 
decrease as l/L^. The constants Be depend 
upon the intrinsic geometry of the loop C but 
not® upon the scale L. 

We have measured the ratios (10) for the 
set of three loops shown in Table 2. Hundred 
measurements have been done every other ten 
of the last thousand iterations of each run of 
the Monte-Carlo simulation described in the 
previous section. Since the runs are done on 
an S'* lattice, finite-size effects are negligible 
only for measurements at scales L = I and 
L = 2. This is enough to get a sensible two- 
parameter fit of the constants A'la^ and Be 
for each loop C separately. Consistency can 
be judged by comparing the values of A'la^ 
between all loops. We find a clear signal that 
the order parameter of the U(2) phase transi- 
tion is indeed the string tension produced by 
the central U(l) subgroup of U(2). The string 
tension A'l jumps discontinuously from zero to 
a non-zero value across the critical curve. A 
typical behavior of this string tension in the 

^We ignore possible logarithmic corrections which 
are very difficult to separate out anyhow. 



critical region is shown in Fig. 6. The value 
of A'la^ is zero within statistical errors every- 
where in the cold phase and for all three loops. 
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Figure 6: String tension A'l of the central U(l) 
subgroup of U(2) in the critical region from the 
loop CI. 



4 U(l) effective charge 

Our introductory analysis of the relation be- 
tween the U(2) lattice gauge theory and 
SU(2)xU(l) suggests that the continuum limit 
should be reached along the lines of constant 
U(l) effective charge in the cold phase. There- 
fore we have looked in this phase for the con- 
tour levels of the fitted parameters Be since 
they should describe the lattice potential in 
the determinant representation. However we 
must resort to some model if we want to la- 
bel the contour levels with a loop-independent 
effective charge to a Coulomb field. A pertu- 
bative calculation has proved to fit the data 
quite well [13]. A first-order expansion gives 
for U(l) 

Be(f3) = ^T(C) + 0(gt) (12) 
47r 

where T(C) is a loop-dependent geometrical 
factor produced by the one-photon exchange 
between parallel edges of the loop. The ex- 
pansion parameter gi is to be interpreted as 



CM 



Table 2: Description of the set of non-planar 
loops included in the measurements. ^ 
means a displacement of ±1 in direction fi. 25 

S(C) is the area of the minimal surface en- 
closed by loop C and T(C) a geometrical fac- q g 
tor described in the text. 

0.15 



7 



an effective coupiing which depends upon the 
coupiing parameters [3 on the iattice. The 
agreement between the fitted vaiues of the ef- 

2 

fective charge ai = ^ for different loops C is 
a stringent test of the applicability of the per- 
tubative approach. These contours are plotted 
in Fig. 7. 



signal-to-noise ratio is of course much higher. 
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Figure 8: Contour levels Vci=0.04, 0.05, 0.06 
and 0.07 (solid lines), Vc2=0.068, 0.086, 0.102 
and 0.120 (dashed lines), Vc3=0.053, 0.068, 
0.080 and 0.095 (dash-dotted lines). 



Figure 7: Contour levels ai=0.12, 0.16, 0.20 
and 0.24 of the effective charge of the central 
U(l) subgroup of U(2) in the cold phase de- 
termined from the loops CI (solid lines), C2 
(dashed lines) and C3 (dash-dotted lines). 

We have taken data in the cold phase on a 
coarse rectangular grid 0.7 < &i < 1.5 , 2.2 < 
&2 < 2.9 with spacing A&i = A&2 = 0.1. 
Nethertheless we can see that the contour lev- 
els are qualitatively the same for all loops 
C even if the roughness of the curves in 
Fig. 7 reminds us that the parametrizations 
of Eqs. (11,12) are only approximate. 

It is also possible to make a model- 
independent analysis of the data and to 
study directly the contour levels of the ra- 
tios Dc{L,j3) or, equivalently, Vc{L,j3) = 
— \nDc{L, [3). We refer loosely to Vc as the 
lattice potential though it is not to be con- 
fused with the static potential. In the deter- 
minant representation, since the string tension 
A'l is zero, we have Vc{L) Be- We exhibit 
in Fig. 8 some contour levels of Vc extracted 
from the potential data at scale 1 where the 



We observe that the contours coincide 
pretty well for all loops, once the level values 
are suitably chosen, and that they are close to 
parallel straight lines, in agreement with the 
naive continuum limit. A closer look shows 
that the parallelism is only approximate and 
that the slope of the contour levels tends to de- 
crease when the potential decreases. For the 
lowest level the slope is about -2.7 which is still 
very different from the naive continuum limit. 
But reaching this limit requires taking hi and 
&2 to cxD and the pertubative corrections are 
known to be large on the lattice. 

We have much more data in the critical re- 
gion and we can attempt to map the contour 
levels in the vicinity of the critical curve in- 
side the cold phase. A typical behavior of the 
effective charge ai across the critical region is 
shown in Fig. 9. It can be seen that a linear fit 
can describe the data near the critical point. 

The values of ai along the critical curve 
inside the cold phase are listed in Table 3 
for all loops we have studied. These values 
are within three standard deviations of their 
average = 0.34. The overall consistency 
among the loops and along the critical curve 
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Figure 9: U(l) effective charge ai in the crit- 
icai region from the loop CI and linear fit in 
the cold phase. 



is pretty good and shows that the pertuba- 
tive parametrization is at least a convenient 
way of describing the data. However we must 
keep in mind that the systematic errors due to 
the hysteresis phenomenon and the pertuba- 
tive parametrization are much larger than the 
statistical errors listed in Table 3. 

The same analysis applied to the potential 
data at scale 1 in the critical region yields 
completely similar results. We shall summa- 
rize these results by stating that our findings 
are consistent with the fact that the first-order 
critical curve is a contour level of the central 
U(l) effective charge and, more generally, a 
contour level of the ratios Dc at scale 1. 

Such limiting values might have a natu- 
ral interpretation in the non-pertubative ap- 
proach [17], based on the Nambu-Goto string 
model, which describes Be as a Coulomb-like 
correction induced by the quantum fiuctua- 
tions of the string. This correction is universal 
in the sense that it does not depend upon the 
gauge group and upon the action of the effec- 
tive string theory which describes the dynam- 
ics of the fiux tube in gauge theories. But it 
does depend upon the dimensionality of space- 
time and the geometry of the loop C. The 
parameters Be have been calculated for arbi- 
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Table 3: Values of the U(l) effective charge ai 
along the critical curve from all loops C . 



trary rectangular loops [18] but not for non- 
planar loops. 

All the results of the last two sections have 
been obtained from measurements at scale 1 
and 2 in lattice spacing units. A noteworthy 
consequence is the smallness of statistical er- 
rors and a very good signal. One might wonder 
how it is possible to extract physical quantities 
from so small lattice distances. But when eval- 
uating the results with respect to conventional 
methods, one should take into account the ab- 
sence of extrapolations, which is the heart of 
our technique. 
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5 Renormalization flow 



We can repeat the analysis of the previous sec- 
tions for the one-scale ratios built with the 
trace of Wilson loops in the fundamental rep- 
resentation of U(2): 



nj 



Tc(L,p) 



(13) 



Again we write 
parametrization: 



< Tr W(C, L) > 
< Ti W(L,L) >^/4 

the same dimensional 



exp{-A'2(S'c 
-B'c + ---} 



(14) 



We have measured the ratios (13) for the same 
set of non-planar loops at the same scales 
L = I and L = 2 and fitted the constants 
K'jd"^ and B'q. Again we find that the string 
tension A'2 jumps discontinuously across the 
critical curve. But now its value stays differ- 
ent from zero in both phases. It is natural 
to associate the non-zero string tension in the 
deconfined U(l) phase to the SU(2) factor in 
U(2). Fig. 10 shows two typical behaviors of 
K'jd"^ in the critical region. 

Along the line hi = 0, inside the cold phase, 
the value of K'jd"^ is very small, K2a^ < 0.015, 
which means that the correlation length ^ = 
l/\/K2 is larger than the lattice, ^ > 8a. This 
finite-size effect explains the inconclusiveness 
of early studies with Wilson action [7] about a 
non-abelian string tension. For hi = 0.4, near 
the critical point, we have K'jd"^ ~ 0.09 and 
a clear signal for a non-abelian string tension. 
However the value of A'2a^ decreases rapidly 
along the line hi = 0.4 and at &2 = 2.8 the 
correlation length ^ is again larger than the 
lattice. Therefore a required step, before at- 
tempting the study of the continuum limit of 
the U(2) lattice gauge theory, is to determine 
the scaling window of our 8^* lattice. The most 
direct method is to use the contour levels of 
K'jd"^ ■ We define the scaling window, some- 
what arbitrarily, by demanding that the corre- 
lation length ^ be smaller than half the lattice, 
K'ja^ > 0.06, to avoid finite size effects. Some 
contour levels are displayed in Fig. 11. 

The same observations that we made about 
the rough appearance of the contour levels of 
the U(l) effective charge apply to Fig. 11. We 
expect the contours to become parallel to the 
hi axis for large hi where one should recover 
the values of the SU(2) string tension with 
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Figure 10: String tension A'2 in the critical 
region from the loop CI and linear fit in the 
cold phase. 



Wilson action. Indeed, in the limit 



the U(l) variables (f) are driven up to a gauge 
transformation into the Z(2) subgroup which, 
by the topological structure of the group U(2), 
is also the center of SU(2). Therefore, in this 
limit, the variables (f) can be absorbed in the 
invariant SU(2) measure and the U(2) model 
with action (6) goes over into the usual SU(2) 
theory with Wilson action. The levels have not 
yet entered this asymptotic regime at hi = l.b 
but they are consistent with the SU(2) data. 

On the other hand, towards smaller hi, the 
contour levels hit the critical curve. We note 
from Fig. 10 that a linear fit can describe the 
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Figure 11: Contour levels A'2a^ = 0.06, 0.12 
and 0.20 in the cold phase determined from 
the loop CI. 



data near the critical point in the cold phase. 
The values of K'jd"^ along the critical curve 
inside the cold phase are listed in Table 4 for 
each loop C . 

The values A'|a^ give the intersection points 
of the contour levels with the critical curve. In 
particular the scaling window of the S'* lattice 
begins, along the critical curve, around &2 ~ 
2.8. We get a good signal-to-noise ratio only 
down to &2 ~ 2.3 where we enter the strong 
coupling regime. 

We are now in a position to test whether the 
renormalization flow of the U(2) lattice gauge 
theory with action (6) is consistent with the 
lines of constant U(l) effective charge. We 
may do it with a reasonable accuracy only 
in the vicinity of the critical curve where we 
can use the values in Tables 3 and 4. Fig. 12 
displays the values of A'^a^ along the critical 
curve on a logarithmic scale as a function of 
&2 = Here go is the bare coupling con- 

slant of the SU(2) gauge held in the continuum 
limit. 

If conflnement persists in the continuum 
limit, the weak coupling behavior of A'2a^ 
should follow the prediction of the renormal- 
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Table 4: Values of the string tension A'2a^ 
along the critical curve from each loop C . 
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ization group: 



non-planar ratios are best suited for such stud- 



2 -^"^'2 
^^2« =A^exp 



(/3offo') 



(15) 



where /3o and fi\ are the usual first two coeffi- 
cients of the /3 function for the SU(2) gauge 
theory and the renormalization scale A de- 
pends upon the ratio t = A(0) is the renor- 
maliztion scale of the SU(2) lattice gauge the- 
ory with Wilson action and the dependence 
A(t)/A(0) is in principle calculable in pertu- 
bation theory. The prediction of Eq. (15) is 
shown in Fig. 12 as the standard band repre- 
senting values of the parameter A in the range 



A 



145 ± 25 



(16) 



We observe that the values A'^a^ extracted 
from the loops CI and C2 agree with their 
average within two standard deviations but 
there is a significant discrepancy with the val- 
ues extracted from the loop C3. The sizeable 
deviation of the loop C3, which explains the 
rather large (estimated) error in Eq. (16), is 
a hint that the systematic errors induced by 
the hysteresis are much larger than the statis- 
tical errors listed in Table 4. Nethertheless the 
string tensions extracted from all three loops 
have the same scaling behavior and there is 
an overall qualitative agreement of our data 
with asymtotic scaling in the expected win- 
dow. Achieving quantitative agreement will 
require not only determining a contour level 
very accurately outside the critical region but 
also either very large lattices or some redefini- 
tion of the coupling constant go [19]. 



6 Universality 



Usually testing universality of a lattice gauge 
model with several parameters requires mea- 
suring a few asymptotic quantities such as a 
string tension K or a glueball mass M and 
checking that their ratios stay constant inside 
the scaling window of the model. Such ob- 
servables are defined in terms of long-distance 
correlations, even in the continuum, and are 
very difficult to extract on the lattice. The 
ratio method of Creutz [16] is a technique 
to study scaling and universality with short- 
distance operators on the lattice from which 
singularities can be removed by renormaliza- 
tion of the bare couplings only. The one-scale 



Let Rc(L, [3) be such a ratio where L is the 
scale in lattice spacing units a and [3 denotes 
all the coupling parameters of the model. In- 
side the scaling window of the model the ratio 
R must satisfy the homogeneous renormaliza- 
tion group equations. Then a curve of con- 
stant ratio define some scale A in physical units 
ofe = A-i/2: 



Rc{L,p{L,X)) = c{X) 



(17) 



Taking the derivative of R with respect to L, 
at fixed A defines a Beta-function. In other 
words the renormalization prescription is to 
set the ratio R at the physical scale A to a 
fixed value c for all values of the cut-off a. 

It follows that, in a lattice gauge model 
which has only one bare coupling in the con- 
tinuum limit, the contour levels of the ratios 
Rc{L,j3) must coincide in the scaling domain 
for all loops C , all scales L and all irreducible 
representations of the gauge group. If there 
are several bare couplings in the continuum 
limit, the contours may depend on what rep- 
resentation defines the ratios Rc{L,j3). We 
have already determined in section 4 the con- 
tour levels of the ratios Dc(L, [3) at scale 1 and 
checked that they do not depend on the loop 
C . Fig. 13 displays some contour levels of the 
lattice potential at scale 1 in the fundamental 
representation of U(2): 



Vc{L,p) = -lnTc(L,(3) 



(18) 



We observe again that the contours coincide 
for all loops C. Our statistics are not quite 
sufficient to check quantitatively that the con- 
tour levels stay the same for the ratios Tc(L) 
at scale 2. But Fig. 11 shows at least a qualita- 
tive consistency with the contours of the string 
tension A'2 which fixes the scale A. 

The contour levels of the ratios Tc(L, [3) dif- 
fer from the contours of the ratios Dc{L,j3). 
This difference is the signature of the decou- 
pling of the SU(2) and U(l) gauge fields. We 
have shown that the contours in the determi- 
nant representation of U(2) characterize the 
U(l) subgroup. We can test whether the con- 
tour levels in the trace representation of U(2) 
characterize the SU(2) subgroup by measuring 
the one-scale non-planar ratios in the adjoint 
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Figure 13: Contour levels Vci=0.16, 0.19, 
0.22, 0.26, 0.31 and 0.37 (solid lines), 
Vc2=0.290, 0.352, 0.413, 0.500, 0.610 and 
0.750 (dashed lines), Vc3=0.232, 0.280, 0.330, 
0.400, 0.485 and 0.590 (dash-dotted lines). 



Figure 14: Contour levels Vci=0.405, 0.475, 
0.54, 0.625 0.73 and 0.845 (solid lines), 
Vc2=0.735, 0.88, 1.02, 1.195, 1.425 and 1.705 
(dashed lines), Vc3=0.595, 0.705, 0.82, 0.955, 
1.15 and 1.35 (dash-dotted lines). 



representation of U(2): 



< TiA W(C, L) > 
< TiaW(L,L) >p/^ 



(19) 



Indeed these ratios depend only upon the 
SU(2) subgroup of U(2). If the continuum 
limit of the U(2) lattice gauge theory is the 
SU(2)xU(l) gauge theory, the adjoint con- 
tours must coincide with the trace contours in 
the scaling region. We exhibit in Fig. 14 some 
contour levels of the ajoint lattice potential at 
scale 1, Vc(L,l3) = - In Ac (L,l3). The dis- 
played level values have been chosen so as to 
make the contours in Fig. 14 coincide pretty 
well with the contours in Fig. 13. 

7 Conclusion 

All the results of our Monte-Carlo simulation 
provide evidence that the U(2) lattice gauge 
theory with action (6) does admit a contin- 
uum limit despite the existence of a first-order 
critical line in the phase diagram. Moreover 
all our data are consistent with the fact this 
continuum limit is the symmetric SU(2)xU(l) 



gauge theory. We have characterized the de- 
coupling inside the cold phase of the central 
U(l) subgroup of U(2) by showing the exis- 
tence of two distinct families of contour lev- 
els. We have shown that the curves of the first 
family are labelled by a dimensionless effec- 
tive charge which can be interpreted as mea- 
suring the coupling to a free massless U(l) 
gauge field in the continuum. We have also 
shown that the curves of the second family 
are labelled by a scale measured in units of 
an SU(2) string tension. We have verified uni- 
versality in the pertubative regime. What we 
have not checked completely is universality in 
the non-pertubative regime. This would re- 
quire more statistics for a better determina- 
tion of the U(l) and SU(2) contour levels at 
larger scales and larger lattices for calculat- 
ing the spectrum. Even if the result is highly 
probable this test is in principle needed to 
prove that the continuum limit of the U(2) 
lattice gauge theory is indeed the symmetric 
SU(2)xU(l) gauge theory. 

We have restricted this study to the upper 
right quadrant (&i > 0,&2 > 0) of the pa- 
rameter plane because we expected, from our 
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preliminary analysis, to find the SU(2)xU(l) 
continuum limit in this domain. However the 
structure of the phase diagram of the model is 
much more complex and extends outside the 
first quadrant. The first-order critical curve 
continues beyond the point 5(0.0,3.3) in the 
quadrant (&i < 0,&2 > 0) where hysteresis 
phenomena still worsen and obscure the in- 
terpretation of the data. Maybe using mul- 
ticanonical algorithms [20] will improve the 
situation. Moreover the existence of another 
critical branch can be deduced as follows [10]. 

Let us consider the partition function of the 
U(2) lattice model with action (6): 

Z(&i,&2) = /" £>[?/] exp ^ &i det [/p 
J p 

jTrt/p + h.c. (20) 

We can choose, on a lattice with free bound- 
ary conditions, a set of links C such that each 
plaquette of the lattice contains exactly one 
link in C. The integration measure is invariant 
with respect to left or right multiplication of 
any link by a U(2) group element. It follows in 
particular that the partition function is invari- 
ant under the change of variable Ui e^'^Ui 
for all links / in £, which implies the identity: 

Z{bi,b2) = Z{bie^'tb2e''^) (21) 

Letting cf) = it, the identity (21) implies that 
the phase diagram is symmetric with respect 
to the hi axis and we can restrict the study to 
&2 > 0. Letting cf) = j we have also 

Z(hi,0) = Z(-hi,0) (22) 

which implies the existence of a critical 
point A' at (hi, h^) = (-1.0,0.0) with the 
same properties as the U(l) critical point 
j4(1. 0,0.0). There must be another first-order 
critical branch which springs from A' , with the 
same slope as the first branch at A, and con- 
tinues towards smaller negative hi. 

We have checked that this second critical 
branch does exist. We find that the gaps de- 
crease along this branch from the U(l) value 
towards zero. Therefore the properties of the 
two branchs are markedly different and there 
is a possibility that the second branch may be- 
come second-order at some point. We hope to 
clarify in another report the nature and order 
parameter of this phase transition and, more 



generally, the properties of the phase diagram 
in the whole upper left quadrant. 

The different global structures of the U(2) 
and SU(2)xU(l) gauge theories can be de- 
scribed in terms of monopoles [12]. There can 
be two types of monopoles in the U(2) gauge 
theory associated to the Z(2) intersection of 
the U(l) and SU(2) subgroups. There has 
been recently much effort [21] devoted to un- 
derstanding the role of monopoles and bound- 
ary conditions in the phase transition of com- 
pact U(l) lattice gauge theory. Enlarging the 
study to U(2) lattice gauge theory would un- 
doubtedly shed a new light on these questions. 
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